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We study the hadron-quark phase transition in the interior of neutron stars. For the hadronic sector, we 
use a microscopic equation of state involving nucleons and hyperons derived within the Brueckner-Hartree- 
Fock many-body theory with reahstic two-body and three-body forces. For the description of quark matter, we 
employ the Dyson-Schwinger approach and compare with the MIT bag model. We calculate the structure of 
neutron star interiors comprising both phases and find that with the Dyson-Schwinger model, the hadron-quark 
phase transition takes place only when hyperons are excluded, and that a two-solar-mass hybrid star is possible 
only if the nucleonic equation of state is stiff enough. 
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PACS numbers: 26.60.Kp, 12.39.-x, 12.39.Ba 



I. INTRODUCTION 

The possible appearance of quark matter (QM) in the inte- 
rior of massive neutron stars (NSs) is one of the main issues 
in the physics of these compact objects. Calculations of NS 
structure, based on a microscopic nucleonic equation of state 
(EOS), indicate that for the heaviest NS, close to the maxi- 
mum mass (about two solar masses), the central particle den- 
sity reaches values larger than 1/fm^. In this density range, 
the nucleon cores (dimension « 0.5 fm) start to touch each 
other, and it is hard to imagine that only nucleonic degrees of 
freedom can play a role. On the contrary, it can be expected 
that even before reaching these density values, the nucleons 
start to lose their identity, and quark degrees of freedom are 
excited at a macroscopic level. 

Unfortunately, it is not straightforward to predict the rel- 
evance of quark degrees of freedom in the interior of NSs 
for the different physical observables, like cooling evolution, 
glitch characteristics, neutrino emissivity, and so on. In fact, 
the other NS components can mask the effects coming directly 
from QM. In some cases the properties of quark and nucleonic 
matter are not very different, and a clear observational signal 
of the presence of the deconfined phase inside a NS is indeed 
hard to find. 

The value of the maximum mass of a NS is one quantity 
that is sensitive to the presence of QM. If the QM EOS is 
sufficiently soft, the quark component is expected to appear 
in NSs and to affect appreciably the maximum mass value. 
In fact the recent claim of discovery of a two-solar-mass NS 
fl] has stimulated the interest in this issue. Purely nucleonic 
EOS are able to accommodate masses comparable with this 
large value ||3-la]. However, the appearance of hyperons in 
beta-stable matter could strongly reduce the maximum mass 
that can be reached by a baryonic EOS 16-101. In this case 
the presence of nonbaryonic, i.e, "quark" matter would be a 
possible manner to stiffen the EOS and reach large NS masses. 
Heavy NSs thus would be hybrid quark stars. In this paper, we 
will discuss this issue in detail. 

Unfortunately, while the microscopic theory of the nucle- 
onic EOS has reached a high degree of sophistication, the QM 
EOS is poorly known at zero temperature and at the high bary- 
onic density appropriate for NSs. One has, therefore, to rely 



on models of QM, which contain a high degree of arbitrari- 
ness. At present, the best one can do is to compare the predic- 
tions of different quark models and to estimate the uncertainty 
of the results for the NS matter as well as for the NS structure 
and mass. Continuing a set of previous investigations using 
different quark models jlll - USIl . we employ in this paper the 
Dyson-Schwinger model (DSM) for QM llTel - fTgtl in combi- 
nation with a definite baryonic EOS, which has been devel- 
oped within the Brueckner-Hartree-Fock (BHF) many-body 
approach of nuclear matter, comprising nucleons and also hy- 
perons. Confrontation with previous calculations shall also be 
discussed. 

The paper is organized as follows. In Sec. |II] we review 
the determination of the baryonic EOS in the BHF approach. 
Sec. [in] concerns the QM EOS according to the DSM, com- 
paring also with the MIT bag model for reference. In Sec. HV] 
we present the results regarding NS structure, combining the 
baryonic and QM EOS for beta-stable nuclear matter Section 
|V] contains our conclusions. 



II. EOS OF HYPERNUCLEAR MATTER WITHIN 
BRUECKNER THEORY 



The Brueckner-Bethe-Goldstone theory is based on a linked 
cluster expansion of the energy per nucleon of nuclear matter 
(see Ref. |20], chapter 1 and references therein). The fun- 
damental quantity of interest in this many-body approach is 
the Brueckner reaction matrix G, which is the solution of the 
Bethe-Goldstone equation, written in operatorial form as 



Gabm=V,b + Y.Y.^ac\pp') 
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where the indices a,b,c indicate pairs of baryons and the Pauli 
operator Qc and energy Ec characterize the propagation of in- 
termediate baryon pairs. The pair energy in a given channel 

c=(BiB2)is 
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7i, (^B, ) + Tb, (kB, ) + Ub, (fefi, ) + Ub2 {kB2 ) (2) 



with 7b (^) = mB + k^/2mB, where the various single-particle 
potentials are given by 
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and are determined self-consistently from the G matrices, 
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The coupled Eq. ([T]i - (|4|i define the BHF scheme with the 
continuous choice of the single-particle energies. In contrast 
to the standard purely nucleonic calculation, the additional 
coupled-channel structure due to hyperons renders the calcu- 
lations quite time consuming. 

Once the different single-particle potentials are known, the 
total nonrelativistic baryonic energy density, £, can be evalu- 
ated: 



=n,p,AX 



k<k 



(B) L 



TB{k) + hjB{k) 



(5) 



It has been shown that the nuclear EOS can be calculated with 
good accuracy in the Brueckner two hole-line approximation 
with the continuous choice for the single-particle potential, 
and that the results in this scheme are quite close to the cal- 
culations which include also the three hole-line contribution 

il. 

The basic input quantities in the Bethe-Goldstone equation 
are the nucleon-nucleon (NN), nucleon-hyperon (NY), and 
hyperon-hyperon (YY) two-body potentials V . The inclusion 
of nuclear three-body forces (TBFs) is crucial in order to re- 
produce the correct saturation point of symmetric nuclear mat- 
ter. The present theoretical status of microscopically derived 
TBFs is quite rudimentary, and in most approaches semiphe- 
nomenological TBFs are used that involve several free param- 
eters usually fitted to the relevant data. An important con- 
straint is the consistency with a given two-body force, i.e., 
both two-body and three-body forces should be based on the 
same theoretical footing and use the same microscopical pa- 
rameters in their construction. Recent results |6, 22] have 
been published within this framework, using meson-exchange 
TBFs that employ the same meson-exchange parameters as 
the underlying NN potential. 

In this paper, we use results obtained in this manner based 
on the Argonne Vjg (V18) [22], the Bonn B (BOB) t24J, and 
the Nijmegen 93 (N93) ll25ll potentials, and compare also with 
the widely used phenomenological Urbana-type (UIX) TBFs 
02611 (in combination with the V\% potential). We remind the 
reader that in our approach the TBF is reduced to a density- 
dependent two-body force by averaging over the position of 
the third particle, assuming that the probability of having two 
particles at a given distance is given by the two-body correla- 
tion function determined self-consistently. 

In the past years, the BHF approach has been extended with 
the inclusion of hyperons |l7|, |27l |28t], which may appear at 
sufficiently large baryon density in the inner part of a NS, 



and lower the ground state energy of the dense nuclear matter 
phase. In our work, we use the Nijmegen soft-core NSC89 
NY potential ll29ll that is well adapted to the available experi- 
mental NY scattering data and also compatible with A hyper- 
nuclear levels 130, 311 . Unfortunately, up to date no YY scat- 
tering data and therefore no reliable YY potentials are avail- 
able. We therefore neglect these interactions in our calcula- 
tions, which is supposedly justified, as long as the hyperonic 
partial densities remain hmited. 

We have previously found rather low hyperon onset den- 
sities of about 2 to 3 times normal nuclear matter density 
for the appearance of the E^ and A hyperons |0, Ufl 123, |2^ 
(other hyperons do not appear in the matter). Moreover, an 
almost equal percentage of nucleons and hyperons are present 
in the stellar core at high densities. The inclusion of hy- 
perons produces an EOS which turns out to be much softer 
than the purely nucleonic case, with dramatic consequences 
for the structure of the NS (see below). We do not expect 
substantial changes when introducing refinements of the the- 
oretical framework, such as hyperon-hyperon potentials ll28ll . 
relativistic corrections, etc.. Three-body forces involving hy- 
perons could produce a substantial stiffening of the baryonic 
EOS. Unfortunately they are essentially unknown, but can be 
expected to be weaker than in the nonstrange sector Another 
possibility that is able to produce larger maximum masses, is 
the appearance of a transition to QM inside the star This will 
be discussed in the next sections. 



III. QUARK PHASE 

The properties of cold nuclear matter at large densities, i.e., 
its EOS and the location of the phase transition to deconfined 
QM, remain poorly known. The difficulty in performing first- 
principle calculations in such systems can be traced back to 
the complicated nonlinear and nonperturbative nature of quan- 
tum chromodynamics (QCD). Therefore one can presently 
only resort to more or less phenomenological models for de- 
scribing QM, and in this paper we illustrate results obtained 
by adopting the DSM. A brief comparison with results from 
the MIT bag model will also be made. 



A. Dyson-Schwinger equations approach 

For the deconfined quark phase, we adopt a model based 
on the Dyson-Schwinger equations of QCD, which provides a 
continuum approach to QCD that can simultaneously address 
both confinement and dynamical chiral symmetry breaking 
ifmllTll . It has been applied with success to hadron physics 
in vacuum llSl [l9l l32l 13311. a nd to QCD at nonzero chemical 
potential and temperature 1134l - [38ll . Recently efforts have been 
made to calculate the EOS for cold quark matter and compact 
stars iMIMl- 

Our starting point is QCD's gap equation for the quark 
propagator 5'(/7;/x) at finite quark chemical potential /i, which 



reads 
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with the renormahzed self-energy expressed as 
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where / represents a translationally invariant regulariza- 
tion of the integral, with A the regularization mass-scale. 
Here, g{pL) is the coupling strength, Dpa{k;ii) is the dressed 
gluon propagator, and Y%{q,p;pL) the dressed quark-gluon 
vertex. Moreover, X" are the Gell-Mann matrices, and m^ is 
the A-dependent current-quark bare mass. The quark-gluon 
vertex and quark wave function renormalization constants, 
Zi 2(C^5^^)^ depend on the renormalization point ^, the reg- 
ularization mass-scale A, and the gauge parameter. 

At finite chemical potential, the quark propagator can as- 
sume a general form with rotational covariance 



S{p;pi 
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i'jpA{p ,pu,u )+B{p ,pu,u 
+ ijAipA + 'M) C{p^,pu, u^) , 



(8) 



where we have written u = (0, ijj.). Please note that we ignore 
quark Cooper pairing herein. Diquark condensate and color 
superconductivity have been considered in the DSM ll4ll - l43ll . 
and we will extend our analysis to that case in the future. 

The kernel, Eq. (jS), depends on the gluon propagator and 
the quark-gluon vertex at finite chemical potential. However, 
little is known about them except at very high chemical po- 
tential, where perturbation theory is applied. We have to ex- 
tend to finite ji the DSM that has been successfully applied to 
hadron physics at /i = 0. 

The Ansatze at zero chemical potential are typically imple- 
mented by writing 



Z,g'Dp^{p-q)Yl{q,p) 



X" 



= mp^qr)DT^ip-q)^^o{q,p) , (9) 



wherein D^p^ {k = p ~ q) 



knka 



{Spa — -p^)A is the Landau- 



gauge free gluon propagator, ^{k^) is a model effective in- 
teraction, and ra{q,p) is a vertex An^afz. 

Herein, we consider the widely used "rainbow approxima- 
tion" 



^a{q,p)^7a, 

and a Gaussian-type effective interaction |14; 
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' In our Euclidean metiic: {■yp,ya} = 2Spa; 7p = Yp', 75 = lAyxfl'n; ah = 
Zt=i Oibr, ab = Y,f=i ajbr, and Pp timelike ^P^ <0. 



involving two parameters D and co. This is a finite-width 
representation of the Munczek-Nemirovsky model Il45i] used 
in Ref. JMXl . which expresses the long-range behavior of 
the renormalization-group-improved effective interaction in 
Refs. ||l9, 46, 47]. Equation (fTTl i delivers an ultraviolet-finite 
model gap equation. Hence, the regularization mass-scale A 
can be removed to infinity and the renormalization constants 
Zi 2 set equal to one. In this model, there is no interaction be- 
tween different flavors of quarks. Therefore, the gap equations 
for the different flavors are independent of each other Here 
we consider the light flavors u, d, and s, neglecting heavier 
ones. 

Usually, there exist two solutions of Eq. ^ in the chiral 
limit, i.e., when nig — 0. One solution with nonzero quark 
mass function M(/?) =B{p)/A{p) is called Nambu solution, 
and represents a phase with dynamical chiral symmetry break- 
ing and confinement. The other solution with zero mass func- 
tion at the chiral limit is called Wigner solution, which repre- 
sents a phase with chiral symmetry and deconfinement. The 
Nambu phase is realized in vacuum, and provides the base- 
ment for describing physics in vacuum. The phase transition 
to deconfinement at finite chemical potential, without consid- 
ering hadron degrees of freedom, is investigated in Ref. ISTII . 
For strange quarks, only the Nambu solution exists in vacuum 
mWi . The appearance of the Wigner phase for strange quarks 
at finite chemical potential will be investigated in the follow- 
ing. 

In Ref. 0371] . the jiq ^ Ansatz is specified as the same as 
in vacuum. It is reasonable in the Nambu phase, but not in the 
Wigner phase. In this article, we introduce a further parameter 
in order to study a possible density dependence of the effective 
interaction. Considering asymptotic freedom at high chemical 
potential, we extend the effective interaction at finite chemical 
potential in a simple manner. 
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introducing the parameter a, which controls the rate of ap- 
proaching asymptotic freedom. In the following we will study 
the dependence of our results on this parameter, marked as 
DSa. Compared to this model, the MIT bag model with non- 
interacting quarks corresponds \a a = °°. For simplicity, we 
assume the same effective interaction Eq. ( fT2] i for each flavor 

The parameters D and ta of the model, Eq. (fT2T i. and the 
quark masses can be determined by fitting meson properties in 
vacuum |44]. We choose the set of parameters O) = 0.5 GeV 
and D=\ GeV^. For simplicity, we use current-quark masses 
niu.d — 0, while OTj = 115 MeV is obtained by fitting the K 
meson mass in vacuum, which is a little different from the 
usually used value mj w 150 MeV in the MIT bag model. 

The EOS of cold QM is given following Refs. L37,^. We 
express the quark number density as 



«-?(]") 
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/«(Ip|;M) = ^y_^a'/'4trDh745^(p;M)], (14) 



where the trace is over spinor indices only. The quark thermo- 
dynamic pressure at zero temperature can be obtained as 



P^i^c) = Pcii^q.a) 
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(15) 



The total density and pressure for the quark phase are given 
by summing contributions from all flavors. For comparison 
with the bag model, we write the pressure of the quark phase 

as 



PQ{llu,lld,lls)= Y^ Pgi^lg) - BdS : 



where 
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(17) 
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Theoretically, we can choose arbitrary values of ^^q, where 
the Wigner phase exists. For massless quarks, the Wigner 
phase exists in vacuum and at finite chemical potentials; the 
results of nq{pLq) are shown in the upper panel of Fig. [1] 
Therefore, we choose /x„.o = Mf/,o = 0, and the correspond- 
ing results of Pqifiq) are shown in the lower panel of Fig. [T] 
For strange quarks, the Wigner phase with nonzero value of 
«s can only be obtained with chemical potential fis and the pa- 
rameter a above some thresholds |48]; see the upper panel of 
Fig. 12] Therefore, we set jx^fl as the value of the starting point 
of the Wigner phase with each a. The single-quark number 
density and pressure P^ for strange quarks are shown in Fig.|2] 
Now only Bqs needs to be fixed in this model. For u and d 
quarks, we can use the "steepest-descent" approximation 114911 . 



P[5]=TrLn[5"^] --Tr[E5] 



(19) 



which is consistent with the gap equation within the "rainbow" 
approximation. In vacuum we obtain the pressure difference 
between the Nambu phase and the Wigner phase for massless 
quarks OSTI] 



; P[Sn] -P[Sw] = 45 MeVfm" 



(20) 



Interpreting the Nambu phase as the real vacuum with P[Sn] = 
0, we then obtain the pressure of the Wigner phase for light 
quarks in vacuum Pi,j{hq = 0) = —45 MeVfm^^ and the ef- 
fective bag constant Bjg — 90 MeVfm^^. However, due to 
the introduction of the parameter a, the 'steepest-descent' ap- 
proximation is not consistent with the gap equation and we 
cannot use it to calculate Pv(/Xj o)- Therefore, in our model we 
simply set Bqs = 90 MeV fm^^ as a parameter, neglecting the 
ambiguity from Ps{jj.s.o)- 

In summary, with respect to the MIT bag model with free 
quarks, the pressure of the DSM at a fixed chemical potential 
is lower, which corresponds to an increasing bag constant as 
chemical potential (density) increases. In the next section, we 
will show a density-dependent effective bag constant for beta- 
stable QM in the interior of NSs. 
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FIG. 1: Quark number density (upper panel) and pressure Pq, 
Eq. l |17t . (lower panel) for massless quarks at finite chemical poten- 
tials. Different curves correspond to different values of the parameter 
a = 0,0.5, 1,2,4,00. The curve denoted DS-MN corresponds to the 
results in Ref I40ll . 



B. MIT bag model 

We adopt the MIT bag model assuming massless u and d 
quarks, s quarks with a current mass of nis = 150 MeV, and 
either a fixed bag constant B = 90 MeVfm^^, or a density- 
dependent bag parameter, 



B(p) = Boo + (Bo - fi») exp [ - j3 (-^' ^ 



(21) 



with Boo = 50 MeVfm-^ Bq = 400 MeVfm^^ and j3 = 0.17. 
This approach has been proposed in lUlll . and it allows the 
symmetric nuclear matter to be in the pure hadronic phase at 
low densities, and in the quark phase at large densities, while 
the transition density is taken as a parameter. Sever al p ossi- 
ble choices of the parameters have been explored in Oil]] , and 
all give a NS maximum mass in a relatively narrow interval, 
1.4M0 <Mmax< 1.7Me. 

It has also been found iH H^l that within flie MIT 
bag model (without color superconductivity) with a density- 
independent bag constant B, the maximum mass of a NS can- 
not exceed a value of about 1.6 solar masses. Indeed, the 
maximum mass increases as the value of B decreases, but too 
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FIG. 2: Same as Fig.[T] but for strange quarks with m^ = 1 15 MeV. 



small values of B are incompatible with a transition density 
p > (2, . . . , 3)po in symmetric nuclear matter, as indicated by 
heavy ion collision phenomenology. (These baryon densities 
are usually reached in numerical simulations |51] of heavy 
ion collisions at intermediate energies without yielding indi- 
cations of "exotic" physics.) For a more extensive discussion 
of the MIT bag model, the reader is referred to I fill . 



IV. RESULTS AND DISCUSSION 

A. Beta-stable hadronic matter 

In order to study the structure of NSs, we have to calculate 
the composition and the EOS of cold, neutrino-free, catalyzed 
matter. We require that the NS contains charge-neutral mat- 
ter consisting of baryons (n, p. A, Yr) and leptons (e^, \i^) 
in beta-equilibrium, and compute the EOS in the following 
standard way \A I SL 15211 : The Brueckner calculation yields the 
energy density of baryon/lepton matter as a function of the 
different partial densities, £{p„,Pp,PA,Pz,Pe,P^i), by adding 
the contribution of noninteracting leptons to Eq. (|5]l. The var- 
ious chemical potentials (of the species / — n,/?,A,E^,e,/x) 
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FIG. 3: Pressure vs the baryon number density of hadronic NS mat- 
ter. Thick curves show results for purely nucleonic matter, whereas 
thin curves include hyperons. 



can then be computed straightforwardly, 

de 
api 

and the equations for beta-equihbrium. 



(22) 



(23) 



{bi and qi denoting baryon number and charge of species /) 
and charge neutrality. 



EP'^'^0' 



(24) 



allow one to determine the equilibrium composition {pi{p)} 
at given baryon density p and finally the EOS, 



' -'-- -^ dp 



e = piJ.„-£. 



(25) 



dp p 

In Fig. [3] we compare the EOS obtained in the BHF 
framework when only nucleons and leptons are present (thick 
lines), and the corresponding ones with hyperons included 
(thin lines). Calculations have been performed with different 
choices of the NN potentials, i.e., the Bonn B, the Argonne 
VI 8, and the Nijmegen N93, all supplemented by a compat- 
ible microscopic TBF |6]. For completeness, we also show 
results obtained with the Argonne V18 potential together with 
the phenomenological Urbana IX as TBFs. 

We notice a strong dependence on both the chosen NN po- 
tential, and on the adopted TBF, the microscopic ones being 
more repulsive than the phenomenological force. The pres- 
ence of hyperons decreases strongly the pressure, and the re- 
sulting EOS turns out to be almost independent of the adopted 
NN potential, due to the interplay between the stiffness of the 
nucleonic EOS and the threshold density of hyperons |6]. The 
softening of the EOS has serious consequences for the struc- 
ture of NSs, leading to a maximum mass of less than 1 .4 solar 
masses laHllfllj which is below observed pulsar masses 1.53(1 . 
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FIG. 4: Baryon density (upper panel) and pressure (lower panel) vs 
the baryon chemical potential of NS matter for different models. The 
vertical dotted lines indicate the positions of the phase transitions 
under the Maxwell construction. 



B. Quark matter in beta-equilibrium 

In compact stars, we need to consider matter in beta- 
equilibrium, 

d^u + e + Ve^u + n + v^^s, (26) 

and charge neutrality, which gives the constraints in pure QM 



fid = llu +lie = liu + Mm = M.V 

2n„ - «</ - «, ^ 



(27) 
(28) 



The numerical results of baryon number density ps = 
{rtu + nii + ns)/3 and pressure versus baryon chemical poten- 
tial Hb = IJ-U + '^■IJ.d are shown in Fig. |4] for the several cases 
discussed above. 

For the hadronic case, we display results only for the Bonn- 
B NN potential, which gives the stiffest EOS without hy- 
perons, and thus is the most favored for reaching large NS 
masses. In this case the thick solid (dashed) curve indicates 



results obtained without (with) hyperons. For QM, we plot re- 
sults with the MIT bag model, with mi,j = 0, m^ =150 MeV, 
and a bag constants = 90 MeVfm^^, or a density-dependent 
B{p), Eq. ( |2TI ) (thin curves). The remaining curves are re- 
sults obtained with the DSM and several choices of the model 
parameter a. 

The crossing points of the baryon and quark pressure 
curves (marked with a square) represent the transition between 
baryon and QM phases under the Maxwell construction. The 
projections of these points (dotted lines) on the baryon and 
quark density curves in the upper panel indicate the corre- 
sponding transition densities from low-density baryonic mat- 
ter to high-density QM. Some qualitative considerations can 
be done. In particular, we notice that the phase transition from 
hadronic to QM occurs at low values of the baryon chemi- 
cal potential when the MIT bag model is used to describe the 
quark phase, whereas much higher values are required with 
the DSM. In some extreme cases, such as DSO, no phase tran- 
sition at all is possible. In fact, the DSM EOS is generally 
stiffer than the hadronic one, and the value of the transition 
density is high. We also notice that with the DSM no phase 
transition exists if the hadronic phase contains hyperons. 

After these indications, we study in the following the phase 
transition with the more sophisticated Gibbs construction. 



C. Phase transition in beta-stable matter 

A realistic model of the phase transition between bary- 
onic and quark phase inside the star is the Gibbs construc- 
tion Jj, [3, ISJIj which determines a range of baryon densi- 
ties where both phases coexist, yielding an EOS containing a 
pure hadronic phase, a mixed phase, and a pure QM region. 
The crucial point of the Gibbs construction, as suggested by 
Glendenning p?], is that both the hadron and the quark phase 
are separately charged, while preserving the total charge neu- 
trality. This implies that NS matter can be treated as a two- 
component system, and therefore can be parametrized by two 
chemical potentials. Usually one chooses the pair (/Xe.,/i„), 
i.e., electron and baryon chemical potential. The pressure is 
the same in the two phases to ensure mechanical stability, 
while the chemical potentials of the different species are re- 
lated to each other, satisfying chemical and beta stability. As 
a consequence, the pressure turns out to be a monotonically 
increasing function of the density. We note that our Gibbs 
treatment is the zero surface tension limit of the calculations 
including finite-size effects |14, 55]. 

The Gibbs conditions for chemical and mechanical equilib- 
rium between both phases read 

Hu + fJ-e = M./ = K«, (29) 

lip + He = M« = Ma = Mm + 2Airf , (30) 

Me-+Mp = 2m«, (31) 

PHilldJ-n) = PQ{lle,lJ-n)=PM{lJ-n)- (32) 

From these equations one can calculate the equilibrium chem- 
ical potentials of the mixed phase corresponding to the inter- 
section of the two surfaces representing the hadron and the 
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FIG. 5: Pressure vs baryon density of NS matter with the Gibbs (up- 
per panel) or the Maxwell (lower panel) phase transition construction 
for different models. 



quark phase, which allows one to calculate the charge densi- 
ties p^ and pc and therefore the volume fraction x occupied 
by QM in the mixed phase, i.e.. 



XP? 



9+{i-x)p!' = o. 



(33) 



From this, the energy density em and the baryon density Pm 
of the mixed phase can be determined as 



Pm = XPQ + {^-X)pH- 



(34) 
(35) 



In Fig. |5] (upper panel) we display results for the EOS in- 
cluding the Gibbs hadron-quark phase transition, using the 
same conventions as in Fig.|4] We notice that the phase transi- 
tion constructed with the DSM turns out to be quite different 
from the one obtained using the MIT bag model. In the for- 
mer case, if the coexistence region does exist, it is shifted to 
higher baryonic density. For comparison, the results with the 
simple Maxwell construction are shown in the lower panel of 
the figure, the main difference being the presence of a plateau 
typical of the first-order phase transitions with one conserved 
charge, and the absence of a mixed phase. 

In order to clarify the fundamental difference between MIT 
model on one side and DSM on the other side, we plot in Fig.|6] 
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FIG. 6: (Color online) Effective bag constant obtained with the MIT 
model and the DSM. 



the effective density-dependent bag constant 



(36) 



obtained with the DSM in NS matter. One observes that B(p) 
is a monotonically increasing function in contrast to the em- 
pirical density dependence introduced with the MIT model. 
We point out that monotonically increasing bag parameters 
have also been found with the color dielectric model lll3[l and 
the Nambu-Jona-Lasinio (NJL) model |56], whereas the den- 
sity dependence in the bag model was introduced by hand in 
order to delay the occurrence of the phase transition as much 
as possible, see Ref. 1. 1 1.1 . It thus appears that an increasing 
bag parameter is a generic feature of many microscopic quark 
models. 



D. Neutron star structure 

We assume that a NS is a spherically symmetric distribu- 
tion of mass in hydrostatic equilibrium. The equilibrium con- 
figurations are obtained by solving the Tolman-Oppenheimer- 
Volkoff equations Il52ll for the pressure P and the enclosed 

mass m. 



dP 

dm 
dr 



Gme (1 +P/£) (1 +Anr'P/m) 



1 — 2Gm/r 



Anre , 



(37) 
(38) 



being G the gravitational constant. Starting with a central 
mass density e{r — Q) = e^, we integrate out until the density 
on the surface equals the one of iron. This gives the stellar 
radius R and the gravitational mass is then 



rR 

MG = m(R)^47t drr^eir). 
Jo 



(39) 



8 



2.5 
2.0 
1.5 
1.0 
0.5 



0.0 



' 1 ' 1 ' 1 
^^^ Maxwell 


1 ' 1 ' 




^^ "^ " " ' 






1 Sibbs 




■ 


/ . - - - - : ; . 


." .\ 


_ 


/' 








/ - .. 


^ - 






/ '' »--'"'"'" 


^rS N 




- 


/ .*. ^rT 


1- \ "^ i. 






V '^ 


n \ ^ 








1 \ \ 


" 


f ,' ! BOB(N+Y) 


: \ t. 




/ ,'•' : DS0.5 


". *'' 




- /..-;/: : - - - DS2 


■■■■■^ -"■■■--•.-.., 


_ 


/ / ...." DS4 


•^:.:.-„V. 




■ t--' MIT90 




•*^i 


MIT-B 

' 1 ' 1 ' 1 


1 ' 1 ' 





0.0 0.6 1.2 1.8 



9 12 

R[km] 



15 



FIG. 7: Gravitational NS mass vs the radius (riglit panel) and the 
central baryon density (left panel) for different EOS employing the 
BOB hadronic model. The thick (thin) lines represent the configura- 
tions calculated with the Gibbs (Maxwell) construction. 



2.0 



1.5 



1.0- 



0.5 



0.0 




Maxwell 
Gitbs 



UIX(N) 

UIX(N+Y) 

DS2 

DS4 



0.0 0.5 1.0 1.5 2.0 10 11 12 13 

P, [fm"'] R [km] 



FIG. 8: Same as Fig.|7l but with the UIX model. 



We have used as input the EOS 
Maxwell construction discussed 
For the description of the NS 
hadronic EOS with the ones by 
in the medium-density regime, 
Metropolis-Teller Il58ll and Baym 
the outer crust. 

The results are plotted in Fij 
play the gravitational mass Mq 
Mq — 2 X lO^^g) as a function 



obtained with the Gibbs or 
above and shown in Fig. |5] 
crust, we have joined the 
Negele and Vautherin lISTIl 
and the ones by Feynman- 
Pethick-Sutherland IH for 

fs. |2] and [8] where we dis- 
(in units of the solar mass 
of the radius R and central 



baryon density p^ . We present results obtained with two ex- 
treme choices of the hadronic EOS yielding rather low or very 
high maximum NS masses, namely, the UIX and BOB mod- 
els, respectively. The maximum masses are in these cases 
1.84 Mq and 2.50 M© with only nucleons, and 130 Mq and 
1 .37 Mq including hyperons. 

The possible effects of the hadron-quark phase transition 
are very different with the MIT model and the DSM: In the 
case of the MIT model, the phase transition begins at very low 
baryon density and thus effectively impedes the appearance of 
hyperons ITJIl . Consequently, the resulting maximum mass of 
the MIT hybrid star is 1.5 Mq, lower than the value of the 
nucleonic star, but higher than that of the hyperon star given 
before. 

On the contrary, with the DSM no phase transition can oc- 
cur and no hybrid star can exist if hyperons are introduced. 
If hyperons are excluded, the phase transition from nucleon 
matter to QM takes place at rather large baryon density. The 
maximum mass of the hybrid star has a slightly smaller value 
than that with only nucleons, and decreases as the density of 
the phase transition decreases. For example, when the nu- 
cleonic Bonn-B NN potential is used, the maximum mass of 
hybrid stars is only a little lower than 2.5 Mq with a = 0.5, 
and decreases to about 2 Mq with a = 2. The same happens 
if the nucleonic UIX interaction is adopted, in which case the 
maximum mass of the hybrid star cannot exceed 1.84M0. We 
remind that with the value a = corresponding to the unmod- 
ified DSM without density-dependent interaction strength, no 
phase transition at all is possible. With a increasing to °o, we 
can obtain a smooth change from the pure hadronic NS to the 
results with the MIT bag model. 

We also notice a large difference in the structure of hybrid 
stars. In fact, whereas stars built with the MIT bag model have 
a pure hadron phase at low density, followed by a mixed phase 
and a pure quark core at higher density, the ones obtained with 
the DSM contain only a hadron phase and a mixed phase, and 
probably no pure quark interior. The scenario resembles the 
one obtained within the NJL model Ill2l[56ll . where at most a 
mixed phase is present, without a pure quark phase. 

Another resemblance to the NJL model is the fact that with 
the Maxwell construction no stable hybrid stellar configura- 
tions are obtained. This is seen in Figs.|7]and[8] where the stel- 
lar configurations built with a Maxwell construction are rep- 
resented by the thin lines, and are characterized by a plateau. 
The sudden onset of the high-density quark phase renders the 
star unstable, with maximum mass values (of purely hadronic 
stars) which turn out to be slightly larger than the ones calcu- 
lated with the Gibbs construction. 

A clear difference between the two models exists as far as 
the radius is concerned. Hybrid stars built with the DSM are 
characterized by a larger radius and a smaller central density, 
whereas hybrid stars constructed with the MIT bag model are 
more compact. 
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V. CONCLUSIONS 

We have investigated the capability of the DSM for QM to 
provide hybrid NS configurations in combination with a mi- 
croscopic hadronic EOS obtained within the BHF formahsm 
including also the appearance of hyperons. We found that the 
unmodified DSM does not allow a hadron-quark phase tran- 
sition, thus requiring the introduction of an empirical density 
dependence of the quark interaction strength. Even in this 
case, however, the early appearance of hyperons in hadronic 
matter inhibits the phase transition. Only in the fictitious case 
of restricting to pure nucleon matter, a phase transition at large 
baryon density is possible, and a hybrid star with 2 Mq is only 
allowed if the nucleonic EOS is stiff enough to produce a NS 
with 2 Mq. Furthermore, stable hybrid stars are only obtained 



with the Gibbs phase transition construction, but not with the 
Maxwell construction. 

These features of the DSM are very different from MIT- 
type quark models, and we have attributed it to the interaction 
between quarks in the DSM, which can be represented as a 
density-dependent bag constant. We have found a different 
density dependence of the effective bag constant in the DSM 
and the MIT quark model. 
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